Abstract 

Intracellular transport of DNA carriers is a fundamental step of 
gene delivery. We present here a theoretical approach to study generi- 
cally a single virus or DNA particle trafficking in a cell cytoplasm. Cel- 
lular trafficking has been studied experimentally mostly at the macro- 
scopic level, but very little has been done so far at the microscopic 
level. We present here a physical model to account for certain aspects 
J of cellular organization, starting with the observation that a viral par- 

^-^ [ tide trajectory consists of epochs of pure diffusion and epochs of active 

f^ ' transport along microtubules. We define a general degradation rate to 

CN| , describe the limitations of the delivery of plasmid or viral particles to 

^ ' the nucleus imposed by various types of direct and indirect hydrolysis 

Q , activity inside the cytoplasm. Following a homogenization procedure, 

which consists of replacing the switching dynamics by a single steady 
state stochastic description, not only can we study the spatio-temporal 
dynamics of moving objects in the cytosol, but also estimate the prob- 
ability and the mean time to go from the cell membrane to a nuclear 
pore. Computational simulations confirm that our model can be used 
to analyze and interpret viral trajectories and estimate quantitatively 
the success of nuclear delivery. 
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Introduction 

The study of the motion of many particles inside a biological cell is a prob- 
lem with many degrees of freedom and a large parameter space. The latter 
may include the different diffusion constants of the different species, veloci- 
ties along microtubules, their number, the geometry of cell and nucleus, the 
number and sizes of nuclear pores, the various degradation factors, and so 
on. The experimental and numerical exploration of this multi-dimensional 
parameter space is limited perforce to a small part thereof, due to the great 
complexity of the biological cell. A great reduction in complexity is of- 
ten achieved by coarse-graining the complex motion by means of effective 
equations and their explicit analytical solutions, which is the approach we 
adopt in this letter. We are specifically concerned with finding a concise 
description of virus and plasmid trafficking in the cell cytoplasm. 

Early vesicle trafficking studies revealed the complex secretion pathways 
[1], whereas much more recent studies of natural (viruses) [2l|3l[ll and syn- 
thetic (amphiphiles) DNA carriers [5] uncover details of the cellular path- 
ways and the complexity of cellular infection. Viruses invade mammalian 
cells through multistep processes, which begin with the uptake of particles, 
cytoplasmic trafficking, and nuclear import of the DNA. However, cyto- 
plasmic trafficking remains a major obstacle to gene delivery, because the 
cytosolic motion of large DNA molecules is limited by physical and chemical 
barriers of the crowded cytoplasm [5 2j • Whereas molecules smaller than 
500kDa can diffuse, larger cargos such as viruses or non-viral DNA parti- 
cles, require an active transport system p^. Viral infection is much more 



efficient than gene transfer using polymers- or lipids-based vectors, where 
a large amount of endocytozed DNA (typically over 100.000 copies of the 
gene) is required to produce a cellular response, while only a few copies seem 
to be necessary in the case of viruses. 

Two recent studies [H [7] showed that microtubules shape the distribu- 
tion of molecular motors and vesicle trafficking inside the cell cytoplasm 
by means of a combination of experiments and numerical simulations. The 
mean concentration of viral species was analyzed in pjj by means of the 
mass-action law. The mechanism of DNA transport in the cytoplasm, how- 
ever, is still an open question. We propose here a coarse-grained reduced 
description of viral trafficking and compare it to plasmid diffusion. Specif- 
ically, we are interested in the probability pjy and the mean time tat for a 
DNA carrier or a virus to get from the cell membrane to a small nuclear 
pore. The evaluation of these quantities calls for a quantitative approach 
to the description of particle trajectories at an individual level and also, to 
quantify the role of the cell organization and the signaling processes involved 
in viral infection. 

We start with the observations that a viral movement can be described 
as a combination of intermittent switches between pure Brownian diffusion 
and active transport along microtubules (figure [T]), while DNA motion can 
be characterized as pure Brownian. We also account for multiple factors 
involved in degradations, such as hydrolyzation, destruction in lysosomes, 
or any other factors that prevent irreversibly the particle from reaching a 
nuclear pore. This degradation process is modeled as killing with a time- 
independent rate k. We use the overdamped Langevin dynamics with killing 
to describe the viral or DNA motion and use Fokker-Planck-type equations 
to obtain asymptotic approximations of p^ and Ti\f in the limit of large and 
small k. We compute the mean time the first among many independent 
viruses reaches a small nuclear pore. Brownian simulations confirm the 
validity of the analytical analysis. The present approach is a first attempt 
to develop a theoretical tool for the analysis of virus dynamics and, hopefully, 
for the study of trafficking of synthetic vectors, a necessary step toward gene 
delivery. 

Modeling Viral or DNA trajectories We model viral trajectories as a 
collection of pieces, each of which is characterized either as directed move- 
ment along microtubules or pure Brownian motion [21 El H]. In contrast, 
DNA motion in the cytoplasm can be adequately described as pure Brown- 
ian motion [9] . Particles moving inside the cell are reflected at impermeable 
surfaces and are absorbed at nuclear pores. A virus travels on microtubules 
as long as it binds to a motor. The three- or two-dimensional position of a 



particle, X{t), is described by the coarse-grained stochastic dynamics 

{y2Dw for a free particle 
(1) 
V(t) for a bound particle 

where w is a (5-correlated standard white noise and V(t) is a time-dependent 
velocity along a microtubule. The velocity V(t) can be either positive or 
negative, depending on whether a viral particle binds to a dynein or to 
a kinesin motor. However, it is not clear what regulatory mechanisms is 
involved in such a choice |12j . 

Mathematical description of a viral trajectory in the cytoplasm. 
We consider the trafficking of a viral particle from an endosome or the cell 
membrane to a small nuclear pore. The cell cytosol is a bounded spatial 
domain i7, whose boundary dVt is the external membrane dQext and the 
nuclear envelope. Most of the nuclear membrane consists of a reflecting 
boundary dNr, except for small nuclear pores dNa, where a viral particle 
can enter the nucleus. We assume that a viral particle that reaches a pore is 
instantly absorbed, so that this boundary is purely absorbing for trajectories. 
The ratio of the surface areas is assumed small, 

'- \dn\ ^ ^^ 

Homogenization of viral trajectory. To replace the intermittent dynam- 
ics between free diffusion and the drift motion along microtubules, described 
in equation ([T]), we use the precise calibration procedure described in [13]. In 
this homogenization procedure, the motion is described by the overdamped 
limit of the Langevin equation 

dX = b (X) dt + V2D dW, (3) 

where D is the diffusion constant and b(X) represents the steady state drift. 
Because all microtubules starting from the cell surface converge to the cen- 
trosome, a specialized organelle located nearby the cell nucleus (figured]), we 
choose in the first approximation a radially symmetric effective drift b(X) 
converging to the nucleus. This approximation can be justified by the study 
[3], where viral trajectories move around the nucleus surface. Although 
viruses move bidirectionally on microtubules, the overall movement is di- 
rected toward the nucleus, thus we only consider here this average motion 
[12| . The homogenized drift in ^ becomes 

b = -i?^ (4) 

where i? is a constant amplitude, which depends on many parameters, such 
as the density of microtubules, the binding and unbinding rates and the 
averaged velocity of the directed motion along microtubules |13j . 
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Figure 1: On the left-side of the ideahzed ceh, a real trajectory consists of 
intermittent Brownian and drift epochs, whereas on the right-side, we show 
two simulated homogenized trajectories. In one of them, the virus arrives 
alive to a nuclear pore, while in the other, it is killed inside the cytoplasm. 
The round dots on the nucleus surface represent the nuclear pores. 

Prom trajectory description to the probability and mean arrival 
time. Viral killing, immobilization or rejection out of the cell and naked 
DNA degradation by nucleases, are coarse-grained into a steady state killing 
rate k. To derive asymptotic expressions for the probability Pat, that a DNA 
carrier (single virus or DNA) arrives to a small nuclear pore alive and for 
the mean time tn, we use the approximation ([2]). The asymptotic estimates 
depend on the diffusion constant D, the amplitude of the drift B, and k. 
These computations are based on the small hole theory [U], which describes 
a Brownian particle confined to a bounded domain by a reflecting boundary, 
except for a small absorbing window, through which it escapes. The domain 
Q. contains a spherical nucleus of small radius 5. The survival probability 
density function (SPDF) p(x, t) to find the virus or naked DNA alive inside 
the volume element x -|- dx at time t is given by 

p(x, t)dyL = Pr{X{t) Gx + dx^r'' > t, t" > t\pi}, (5) 

where r'^ is the first passage time of a live DNA carrier to the absorbing 
boundary dNa, r is the first time it is hydrolyzed, and pi is the initial 
distribution. The SPDF p(x, t) of the motion ([3]) is the solution of the mixed 



initial boundary value problem for the Fokker-Planck equation (FPE) [15] 
^(x,t) = L>Ap(x,t)-V-b(x)p(x,t)-Mx,t) 
p(x, 0) = Pj(x) for xGfi 
p(x, t) = for X G dNa 

J(x, t) • Hx = X e dNr U Sr^exi, 

where nx is the unit outer normal at a boundary point x. The flux density 
vector J(x, f) is defined as 

J(x, t) = -DVp{^, t) + b(x)j9(x, t). (6) 

The probability P/v that a live DNA carrier arrives at the nucleus is P/v = 
Pr{T^ < T^} [in]- This probability can be expressed in terms of the SPDF 
[16] by 



Pjv = 1 - Pr{T^ >T^] = 1 - / A;(x)p(x) dx, 

Jo. 

where p(x) = J^ p(x, t) dt is the solution of equation 

qL>Ap(x) - V •b(x)p(x) -fc(x)p(x) = -pi(x) for x G 17 (7) 

q with the boundary conditions ©• Using the pdf of the time to absorption, 
conditioned on the event that the DNA carrier escapes alive Pr{r° < i | r"^ < 
r'^}, we define the conditional mean time to absorption as 

rjv = EIt" It" <t^]= (1 - Pr{r'^ Ktlr" <t'') dt. 

Jo 

Following the computations of [ITj, we get 

^ /q^(x)q?x - /^ A;(x)g(x)dx 

1 - /^ A;(xMx) dx " ^'^ 



where 



satisfies [T7] 



/>oo 

g(x) = / sp(x, s) ds (9) 

JO 



- p = DAq{yL) -V -hqj-kq for x G J7 (10) 

with boundary conditions ([6]). 
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Asymptotic expressions: the plasmid case. To obtain explicit expres- 
sion for P/v and tn for a nucleus containing n well separated small holes 
(nuclear pores) on its surface, we consider first a killing rate k smaller than 
the diffusion constant D. The asymptotic analysis for naked DNA (b = 0) 
leads to [IT] 
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where ^ = tttt /q ^(^) "^^i and rj is the radius of a small absorbing disk. 
Formula (jlip does not depend on the specific shape of the killing rate k, but 
rather on its integral. We compare this asymptotic formula with Brownian 
simulations obtained for parameters R = 20fim; 6 = ^; r] = Sj^ = 1.05/im; 
^ = 3m^~^ [E]; £» = OMfim'^s-^ [9]; n = 1, (a single big hole). This 
simulation corresponds to a cell with 2% of the nuclear surface occupied by 
a large nuclear pore, or equivalently, to a simulation with n = 2000 pores 
of radius 25ra?7i[19]. Because formula ([TT]) depends only on the product rnj, 
both simulations give the same result. 



Time and Probability 


TN 


Pn 


Theoretical values 


3567s 


0.90% 


Simulated values (2000 particles.) 


3564s 


0.97% 



When A;(x) is much larger compared to diffusion, a boundary layer analysis 
leads to the asymptotic expression 
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where we assume that the smooth killing rate is a constant k^ in the neigh- 
borhood of the nuclear surface. Asymptotic expressions: the virus 
case. For a virus trajectory governed by equation ([3]), with a constant 
scalar drift B, the leading order term of the probability and the mean time 
are given by [T7] 



^(l^6^ + 2(^?6 + 2(^) 
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^A§6^ + 2{^Y6 + 2{§r]k + l ^Jp^ + 2{§y6 + 2{§r]k + l 



These asymptotic formulas show that the main contribution to the prob- 
ability and the mean time comes from a boundary layer located near the 
nucleus surface. The killing rate k in this case is the averaged value of 
the killing field in that boundary layer |17j . In figure O we compare the 
probability to arrive alive at the pore and the mean arrival time for several 
values of the drift and the constant killing rate. For a large number of 
microtubules, the drift B equals the apparent velocity [T3] (10% [20] of the 
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Figure 2: MFPT (top left) and the arrival probability (top right) for increasing 
values of the drift with k = j^s"^ and for increasing values of the steady state 
killing rate {B = 0.2/ims~^) (bottom). 2000 trajectories are simulated, theoretical 
values and simulated ones are drawn with dashed and solid lines, respectively. 
R = 20/xm; 5 = Aiim; V = jz^ = l-OS/Ltm; D = 1.3/Lim^s~^[l]; n = 1. 
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Figure 3: MFPT 

of the first virus. 
300 trajectories are 
simulated (solid 
line). R — 20 fim 

S = 4/i7TT. 

rj = j^S = 1.05/1771 



1.3/im^ 



D = 

n 

B = 0.2fims 



1 and 



minus end velocity, approximatively equal to 2^m/s [4J). Using formula 
(|12p . we can now predict the effect of changing the effective drift B = 0.2 



by ±30%: increasing the drift leads to a probability P^ = 0.80 and a 
mean time t^ = 731s, while reducing the drift gives P^ = 0.64 and 
Tj^ ° = 1293s. We conclude that decreasing the drift increases the time 
by 33% (ttv = 974s) and decreases the probability by 12% {Pn = 0.73), 
while increasing the drift, reduces the time by 22% and increases the prob- 
ability by 10%. These results show the nonlinear effect of the drift. In a 
biological context, decreasing the drift can be implemented by disrupting 
the microtubule network. 

Mean first passage time of the first virus to reach the nucleus. 
When M viruses enter a cell, the number Ma of live viruses arriving at the 
nucleus is given by Ma = PnM. The mean time the first live virus arrives 
at a nuclear pore is given by 

\ A:=0 fc=l / 

where xi = 1 — P/v- Finally, asymptotic expansions give 

' 2^ (l + ^xi) for xi « 1 

for xi ~ 1 and M » 1. 
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The theoretical results are compared with Brownian simulations in figure [3l 
The closed form expressions II HfT2] facilitate the exploration of the multi- 
dimensional parameter space of cellular delivery of both DNA and virus 
trafficking. Cytoplasmic trafficking is a limiting step of gene delivery. Eluci- 
dating viral motion in the cytoplasm may provide a quantitative tool for the 
improvement and optimization of delivery of synthetic vectors. The present 
approach can provide a resource for optimizing the design of synthetic vec- 
tors and for the analysis of the parameters of viral infection. 
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